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THE  EFFECTS  OF  TRANSVERSE  GRADIENT  OF  STATIC  MAGNETIC  WIGGLERS 

ON  THE  FREE  ELECTRON  LASER 

I.  Introduction 

The  operation  of  the  free  electron  laser  (FEL)  requires  the  interaction  of  a 

relativistic  electron  beam  with  the  electromagnetic  field  combined  with  the 

wiggler  field.  Most  of  the  FEL  experiments  in  progress  utilize  periodic  static 

magnetic  field,  B^,  which  can  be  either  linearly  polarized  or  circularly 

polarized  (helical  magnetic  field).  The  unrealizable  one-dimensional  model  of 

the  wiggler  field  is  adequate  for  the  study  of  the  fundamental  phenomenons  of  the 

FEL  process.  For  more  refined  calculations,  the  realizable  wiggler  field,  which 

satisfies  V  •  B  -  0  and  V  x  B  -0,  will  have  to  be  considered.  Ihe  realizable 
~w  ~w 

wiggler  has  transverse  variation  as  well  as  a  non-zero  axial  component  of  the 
magnetic  field  off  the  axis. 

In  this  paper,  we  derive  the  approximate  equations  of  motion  in  a  realizable 
circularly  polarized  wiggler  field. **2  Those  particle  trajectories  are  used  to 
obtain  effects  due  to  transverse  gradients.  We  will  compare  the  results  against 
that  due  to  a  realizable  linearly  polarized  wiggler  field  derived  elsewhere"*'^. 
The  axial  component  of  the  magnetic  field  causes  the  betatron  oscillations.  The 
betatron  oscillations  of  an  electron  beam  possessing  finite  emittance  results  in 
an  effective  axial  velocity  spread  in  the  axial  direction  for  an  originally  cold 
electron  beam  possessing  finite  radius.  We  find  that  the  requirement  for  the 
power  density  of  input  radiation  field  to  initially  trapping  all  the  electrons 
with  a  circularly  polarized  wiggler  is  one  quarter  of  the  necessary  input 
radiation  field  for  an  equivalent  linearly  polarized  wiggler. 

Manuacript  lubmitted  April  2,  1982. 


II.  Particle  Trajectories  for  a  Helical  Higgler 

The  circularly  polarized,  static,  magnetic  field  of  a  contoured  wlggler  can 
be  written  as  the  gradient  of  a  scalar  potential 

-2B 

$(r,0,z)  ~  — —  IjC^r)  co8(kwz  -  0),  (1) 

w 

where  B  »  -V4>,  I,  is  the  modified  Bessel  function  of  the  first  order,  B  is  the 
~w  1  w 

amplitude  of  the  magnetic  field  on  axis  and  is  the  wavelength  of  the  magnetic 
field  on  axis,  i.e.,  r  =  0.  The  corresponding  magnetic  field  is 

5W  -  IUIo(kwr)  cos(kwz)  +  I2(kwr)  cos(kwZ  "  20>]  ex  (2) 

+  BwfXo(kwr)  sln<kwz)  _  X2(kwr)  sin(kwZ  "  20>]  ey 
‘  2  Bw  Il(kwr)  sln(kwz  '  0>  V 

For  the  calculation  of  particle  trajectories,  it  is  more  convenient  to  write  the 
magnetic  field  in  terms  of  the  vector  potential  B#  »  V  n  Aw,  where 

Aw  2+  Aw[-X0<kwr)  cos(kwz)  +  I2(kwr)  cos(kwZ  ~  20>]  ex 

“  A  [i  (k  r)  sin(k  z)  +  I,(k  r)  sin(k  z  -  20)]  e 
w  ow  w  z  w  w  y 

and  A  ■  B  /k  . 
w  w  w 

We  will  solve  the  equations  of  motions  by  making  a  perturbation  expansion  of 
the  transverse  momentum, 
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P  «  P  +  6P  , 

XX  X 


P  -  P  +  6?  , 

y  y  y 


where  P  ^  »  6P  and  P  »  6P  .  Furthermore  we  will  assume  the  electron 
x  x  y  y 

beam  to  be  confined  sufficiently  close  to  the  center  of  the  axis, 

i.e.,  1  »  k  r  and  assume  v  »  v  ,  v  .  Since  B  »  B  and  B  »  B  ,  we  will 
w  z  x  y  x  z  y  z 

define 


^  v  B  , 
dt  c  zo  y 


—JL —  -  v  B  ,  (5b) 

dt  c  zo  x* 

Y  2  -  1 

t  ©  *•  2  1/2 

where  v^q  ■  c  [ - - ]  is  the  axial  velocity  for  the  particle  that 

Yo 

g 

remains  on  the  axis  and  Bq^  =  ( — (j-^)  is  the  transverse  velocity.  Vfe  will 

Y  m  c  w 
o  o 

0 

consider  the  steady  state  limit,  where  ■  0.  Making  use  of  the  fact  that 


vzo  If  »  vx  h  +  vy  I?  we  can  wrtte  It  i  vzn  If*  To  the  lowe8t  order*  we  find 


dt  —  zo  3 z 


p  <°>  ^  i±L  (A  (?  e,  Z)  .  e  ), 

v  —  r  v  's.u  ’  9  v  *  9 


-  Z)  *  ey)» 


JL 


where  r,  0  are  functions  of  the  Lagrangian  variables  (xq,  yQ,  vxo,  VyQ»  z)» 
xq,  yQ  are  the  initial  transverse  coordinates  at  z  ■  0,  and  vxq>  v^q  are  the 
initial  transverse  velocities  at  z  ■  0.  In  the  limit  where  r  ■  0,  we  recover  the 
one-dimensional  result. 

The  equations  that  govern  the  next  order  corrections  are: 


d<SP 
_ x 

dt 


B 

yo  z 


-  B  6v 
y  z 


), 


B 


x 


v  B  ), 
xo  z  * 


(7a) 

(7b) 


-(k  r)2 

where  6v  *  - ^ -  v  3  ,  [2-cos  2(k  z  -  0)].  At  this  point,  we  perform  a 

z  8  zo  ol  w 

transformation^-^  from  the  Eulerian  independent  variables  to  the  Lagrangian 

independent  variables  (x  ,  y  ,  v  ,  v  ,  z).  The  total  derivatives  in  time  for 
v  o’  Jo  xo  yo 

the  left-hand-side  of  Eqs.  (7a, b)  can  be  rewritten  as 


d6P 
_ x 

dt 

d6P 


r  d  X  .  z  dxi 
ymv  Iv  — r- +  j — — I , 

o  o  zo  1  zo  ^2  dz  dzJ 

,2~  dv  - 

y  m  v  [v  — ?■  +  ~t~-  , 

o  o  zo1  zo  ^2  dz  dz J 


(8a) 

(8b) 


where 


-  2k  c  B  ,2I,(kJ)  I(kr)  sin  2(k  z  -  0). 
dz  w  ol  2  w  o  w  w 

We  look  for  the  slowly  varying  solutions,  where  ^  «  kwx  and  &  «  k^y 

for  k^r  «  I.  In  the  slowly  varying  limit,  it  is  a  good  approximation  to  assume 


that  the  expressions,  which  contain  sinusoidal  terms  of  0<wz)>  contribute  little 

^  /Si  _  ^ 

to  r  and  6.  The  appropriate  equations  for  x  and  y  become 


(9a) 


2~ 

=  0  ?  k  2  I,(k  r)[l  (k  r)  cos  0  -  I_(k  r)l 
.2  oi  w  1  w  1  o  w  2  w  J 

dz 


=  0  2  k  2  I  (k  r)  I  (k  r)  sin  0. 

,2  oi  w  1  w  o  w 
dz 


Taylor  expanding  the  Bessel  functions  and  keeping  only  the  lowest  order  term, 


Eqs.  (9a, b)  reduce  to 


—  =  K  2x 

a  2  o  ' 

dz 


(10a) 


&  -  K  .  2y. 
,  l  o 
dz 


(10b) 


where  K  *  0  ,k  //2  is  the  wavenumber  of  the  betatron  oscillation.  We  have  used 
o  oi  w 

the  fact  x  =  r  cos  0  and  y  »  r  sin  0.  The  initial  conditions  are  x(z  *  0) 

=■  x  ,  y(z  =  0)  =  y  ,  ~|  =  -22.  and  ^-|  =  -^2. 

o  o  dz  z=o  v  dz  z=o  v 

zo  zo 

We  find  that  electrons  execute  betatron  oscillations  in  realizable  wigglers 
with  transverse  gradients.  The  solutions  to  Eqs.  (10a, b)  are  simply 


cosK  z 
o 


- —  sin  K  z  x 
K  c  o  o 

o 


v  »  I  -K  c  sinK  z  cos  K  z I  I  v  , 
x  o  o  o  xo 


(Ha) 

(lib) 


y  *=  l  cosKqz 


— —  sin  K  z  y 
K  c  o  Jo 

o 


v  “  I  -K  c  sinK  z  cos  K  z [  v  I  . 
yj  L  o  o  o  J  L  yoj 


(12a) 

(12b) 


An  informative  alternate  approach  to  write  Eqs.  (11a)  and  (12a)  is 


6 


JL 


X  «  Xg  COs(KqZ  -  #x)» 

y  -  yB  cos(Kqz  -  *y). 


(13a) 

(13b) 


where  xfi  -  [(-| - ~)2  +  xq2]^2  and  yfi  =  [(^ - ^p)2  +  y^]^2  are  the  amplitudes 


-1  Xo  -1  yo 

of  the  betatron  oscillations,  and  $  =  cos  ( — )  and  ♦  ■  cos  ( — )  are  the 

x  XB  y  yB 

initial  phases  of  the  betatron  oscillations  in  the  x  and  y  direction 
respectively.  The  transverse  velocities  associated  with  the  betatron 
oscillations  are 


v  =  -K  c  x.  sin  (K  z  -  $  ), 
x  o  Ti  o  x  ’ 


v  =  -K  c  y_  sin(K  z  -  0  ). 
y  o  B  o  y 


To  summarize,  the  particle  trajectories  are 


(14a) 


(14b) 


P - -  (A  (r,  0,  z)  •  e  )  +  Y  m  v  , 

x  —  c  '~w  ’  ’  x  o  o  x 


(15a) 


P  ~  - —  (A  (r,  9,  z)»e)+Ymv. 
y  —  c  ~w  y  o  o  y 


(15b) 


III.  Emittance  and  Axial  Velocity  Spread 


All  realizable  electron  beams  have  finite  emittance.  The  emittance  in  the  x- 

direction,  e  ,  is  defined  as  e  *  v  x  (v  /c  )  where  x 

’  x  x  o,raax  ro,max  o,max 

and  v  are  the  electron  beam  waist  and  the  maximum  velocity  spread 

xo,max 

respectively  in  the  x-direction  at  the  entrance  of  the  magnetic  wiggler.  The 

condition  for  the  envelope  of  the  electron  beam  to  remain  approximately  constant 

inside  the  realizable  wiggler  with  transverse  gradient  is  to  require  that 

particle  (x  v  =  0)  and  particle  (x  =  0,  v  )  have  identical 

o.raax,  xo  o  xo,max 

betatron  amplitude,  i.e., 


x  =  v  /(K  c)  =  (e  /(itK  ))1/2. 

o,max  xo,max  o  x  o 


(16) 


The  condition  on  x  in  (16)  also  leads  to  the  smallest  electron  beam  radius 

o.raax 

inside  the  wiggler  field, 


2  e 


r— 


xi  1/2 


Similarly,  we  will  require 


(17) 


o 


(18) 


where  e  =  iry  (v  /c)  is  the  emittance  in  the  y-direction,  y  and 

y  7o,max  yo,max  o,max 

v  are  the  electron  beam  waist  and  the  maximum  velocity  spread  respectively 

yo.max 

in  the  y-direction. 

The  gradient  in  the  wiggler  field  will  lead  to  an  increase  in  the  energy 

spread  associated  with  the  axial  motion.  Let  us  consider  a  cold  electron  beam 

2 

with  total  energy  (y  -  1)  mQc  .  Since  part  of  the  energy  is  associated  with  the 


m 


"i 


transverse  motion,  the  axial  velocity  of  a  particle  decreases  as  the  amplitude  of 
the  betatron  oscillation  increases.  Finite  emittance  in  the  transverse  direction 
contributes  to  effective  energy  spread.  The  axial  velocity  can  be  written  as 


(P  (0)  +  6P  )2  +  (P  (0)  -HSP  )2 

'  y  ^  y  y 

(y  m  c)2 
oo 


We  are  interested  in  the  mean  axial  velocity.  Thus,  we  will  drop  terms  that  are 

sinusoidal  functions  of  k  z.  The  effective  axial  velocity  is  found  to  be 

w 


<v(x,y,v  ,  v  )  >  =  v  -  Av„  (x  ,  y  ,  v  ,  v  ), 
z  o*  o  xo  yo  zo  Boo  xo  yo 


(19) 


Y  2  -  1 

T  *■  2  1  /  2 

where  v  =  c  f - = - 3  .l1  is  the  axial  velocity  in  the  absence  of  the 

zo  1  2  ol  J 


betatron  oscillations  in  a  non-realizable  wiggler  field  without  transverse 

gradient,  and  is  also  the  axial  velocity  of  the  electron  travelling  on  axis  in  a 

realizable  wiggler  field,  <(  )>  m  j-  /L(  )dz  is  an  average  over  axial  distance  L 

Li  O 

for  L  »  and  Av  is  the  change  in  the  axial  velocity  for  a  particle 

k  p 

w 

executing  betatron  oscillation.  Since  the  axial  velocity  spread  is  proportional 
to  the  radial  excursion,  the  maximum  axial  velocity  spread,  due  to  betatron 
oscillation  alone,  is 


Av, 


2c(8olVb/2^ 


(/2/ir)c  BoJkwe, 


(20) 


where  r.  ■  x,  -  y.  is  the  radius  of  the  electron  beam  and  e  -  e  ■  e  is  the 
b  d  b  x  y 

emittance  for  an  axially  symmetric  beam.  The  corresponding  longitudinal  energy 

2  2 

spread  is  related  to  the  velocity  spread  by  AE^  =  yzq  (AVg/c)YomQc  .  The 
longitudinal  energy  spread  in  terms  of  the  radius  of  the  electron  beam  and  the 


8 


emittance  are  spectively 


AE, 


Y  d  ,k  r,  0 
^  zo  ol  w  bj2 


Y  in  c 
o  o 


(/2/n )  (  8  k  E. 

zo  ol  w 


(21) 


One  efficiency  enhancement  approach  is  to  initially  trap  a  large  fraction  of 
the  electrons  in  the  ponderomoti ve  potential  well  and  adiabatically  extract 
kinetic  energy  from  the  particles.  In  order  to  trap  a  substantial  fraction  of 
the  electrons,  we  require  the  trapping  potential  to  be  larger  or  at  least 
comparable  to  the  axial  energy  spread,  i.e.,  |  e  1 4>t  rap  >  AE^.  The  initial  depth 
of  the  trapping  potential  is 


2  1/2 
le|*trap/(YomoC  >  =  4Yoz£?oi(AR/A„)  * 


(22) 


Combining  Eqs.  (20),  (21)  and  (22),  we  obtain  the  requirement  for  the  input 
radiation  field  in  terms  of  the  radius  of  the  electron  beam. 


\  >  (Yzo 


2  2  2 
B  .  k  r.  ) 

ol  w  b 


(23) 


Since  the  optimal  radius  is  determined  by  the  emittance,  we  can  also  rewrite  the 
input  radiation  field,  required  for  trapping  the  electrons,  in  terms  of  emittance 


A  ^  r  zo  w  £-|2  w,Cir 

Vdr  >  - S - i  T  * 


(24) 


Now,  we  will  compare  the  effects  of  transverse  gradient  on  identical  and 
axially  symmetric  electron  beams.  When  the  transverse  gradient  is  not 
considered,  we  have  previously  found  that  the  FEL  action  utilizing  a  linearly 
polarized  wiggler^  and  a  circularly  polarized  wiggler  are  qualitatively 


equivalent  if  A 

w,Cir 

we  notice  that 


Aw,Lin//2  and  Vcir  “  Vun^2* 

Ko,Cir  =  Ko,Lin//2’ 

r  =  21/4  r 

b.Cir  b.Lin’ 


Using  these  relations. 


(25) 


and 


AYB,Cir  “  AY6,Lin//2> 


^trap.Cir  ^trap.Lin^2’ 


The  wavelength  of  the  betatron  oscillation  is  longer;  the  minimum  radius  of 
the  electron  beam  in  the  wiggler  is  larger  and  the  axial  velocity  spread  is 
smaller  for  the  circularly  polarized  wiggler  than  for  the  linearly  polarized 
wiggler.  The  mimitnum  radiation  field  required  for  trapping  the  electrons  for  the 
linearly  polarized  wiggler  is 


,  ^  cYzo  kw  e->2  Aw,Lin 

R,Lin  tt  4  • 


(26) 


The  ratio  of  the  input  radiation  power  densities  required  for  trapping  is 


P 


Cir 


_  p 

4  Lin’ 


(27) 


-12  2  -12  2 

where  PCir  “  (4itc)  u  A^  and  P^^  *  (8itc)  a)  AR  L^n»  The  minimum  initial 
power  density  required  for  trapping  the  electrons  for  a  circularly  polarized 
wiggler  is  only  one  quarter  of  that  for  a  linearly  polarized  wiggler.  Therefore, 
if  the  emittance  combined  with  the  betatron  oscillation  is  the  largest  source  of 
energy  spread,  then  the  circularly  polarized  wiggler  is  a  better  choice  for 
axially  symmetric  electron  beams. 


10 


IV.  Numerical  Results 


We  illustrate  the  effect  of  transverse  gradients  within  a  realizable 

magnetic  wiggler  with  the  following  numerical  example.  We  consider  a  wiggler 

with  a  period  l  =  2.8  cm  (k  =  2.24  cm  ^),  and  magnetic  field  B  =  5  kG 

w  w  w }  Li  n 

3 

(A^  =  2.23  x  10  statvolts)  for  a  linearly  polarized  wiggler  or  Bw 

3 

=  5/2  kG  (A^  =  1.58  x  10  statvolts)  for  a  circularly  polarized  wiggler;  and 

we  consider  an  electron  beam  with  energy  of  25  HeV  (y  =  50)  and  e  =  2it  mm  mr. 

The  parameters  are  appropriate  for  radiation  at  10.6  pm.  The  wavelength  of  the 

-2  -1 

betatron  oscillation  is  76  wiggler  periods  (K  =  2.9  x  10  cm  )  for  the 

o ,  Lir 

_2 

circularly  polarized  wiggler  and  is  54  wiggler  periods  (K  ,  .  =  4.1  x  10 

o.Lin 

cm  S  for  the  linearly  polarized  wiggler.  We  find  that  the  optimal  radii  of  the 

electron  beam  within  the  wiggler  field  are  r,  =  0.12  cm  and  r.  TJ  =0.1  cm. 

b,Lir  b, Lin 

These  satisfy  our  approximation  k^r^  «  1.  The  input  radiation  power  densities 

2  2 
required  for  trapping  are  =  25  MW/cm  and  P  “  97  MW/cm  . 

First,  we  will  examine  the  accuracy  of  the  expression  for  the  particle 

trajectories,  Eqs.  (15a,  b).  We  will  define 


a 

x 


=  [P 

L  x, exact 
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and 
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where  P  _  and  P  _  are  particle  moraentums  in  the  x  and  y-directions 

x, exact  y, exact 

obtained  numerically  using  the  exact  wiggler  magnetic  fields  Eq.  (2),  and  P^q  is 
the  total  initial  transverse  momentum.  Figure  1  shows  plots  of  average 

<N, 

o  and  o  as  a  function  of  axial  distance  z  for  a  particle,  that  is  initially 
x  y 

ii 


•  '  W  ft*.'  ’  0! 
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at  x  ■  .12  cm  and  y  -  v  ■  v  “0.  For  this  example,  k  r.  ~  0.25,  and  the 

error  is  only  1%.  Next,  we  illustrate  betatron  oscillations  for  k^r^  comparable 

to  1.  Figure  2  contains  plots  of  x  as  a  function  of  axial  distance  z  for  various 

values  of  k  r.  with  v  =  v  =  0.  We  find  that  as  k  r.  increases  the  wavelength 
w  b  xo  yo  w  b 

of  the  betatron  oscillation  decreases,  while  the  amplitude  stays  constant.  Thus, 
the  effective  axial  velocity  spread  becomes  larger  than  the  expression  shown  in 
Eq.  (20)  as  k  r,  becomes  larger. 


kwrb=l 


kwrb= 


AXIAL  DISTANCE  Z 

Fig.  2  —  Plots  of  x  as  a  function  of  axial  distance  z  for  particles 
with  yG  =  vxo  =  vvo  =  0  and  various  values  of  kwrb 
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Appendix  I 


Betatron  Oscillations  Due  To  Transverse  Gradient 
Of  The  Linearly  Polarized  Wiggler  Field 

This  appendix  will  obtain  the  particle's  transverse  motion  and  the  axial 
velocity  spread  inside  a  linearly  polarized  realizable  wiggler  field.  For  the 
illustrative  purposes  here,  we  will  assume  that  the  wiggler  amplitude  and 
wavelength  are  constant  in  the  axial  direction.  The  vector  potential  of  the 
wiggler  field  is 


A  (y,z)  «  A  cosh(k  y)  cos(k  z)  e  .  (1.1) 

~w  J  w  '  \r  w  x 

The  associated  wiggler  magnetic  field  has  a  z  component  off  the  axis,  in  addition 
to  the  y  component, 

Bw(y,z)  »  -A^lcoshCk^)  sin(kwz)  ey  (1.2) 

+  sinh(kwy)  cos(kwz)  ez]. 


For  particle  trajectories  in  the  transverse  directions,  it  is  sufficient  to 
assume  that  their  motion  is  effected  by  the  wiggler  field  alone,  i.e.,  Px(y,z) 


The  particle  location  in  the  x-direction  can  be  obtained  by  Integrating  the 
momentum  in  the  x  direction. 


6  .  v 

^  o  X  ^  xo 

x  ~  x  +  k —  co8h(k  y)  sin  k  z  + - z 

°  w  W 


(1.3) 
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where  x  and  y  are  functions  of  (z,  x  ,  y  ,  v  ),  x  and  y  are  the  initial 

o  o*  yo  *  o  o 

transverse  coordinates,  v^o  is  the  initial  transverse  velocity  in  the  y- 
direction. 

The  particle  motion  in  the  y-direction  is  due  to  the  z  component  of  the 
wiggler  field. 


v  B  . 


x  z 


(1.4) 


We  will  assume  that  the  fast  oscillatory  terms,  with  wavelength  half  the  wiggler 
wavelength,  are  unimportant,  and  replace  vg  by  c  at  appropriate  places.  We  find 
that 


H  +  ¥kw slnh  2V  - 0 


(1.5) 


where 
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This  equation  can  be  integrated  once  to  give 


d£ 

dz 


P  v 

(-{-  (cosh(2kwyo)  -  cosh  (2kwy))  +  (~^p-)2) 


1/2 


(1.6) 


To  find  an  explicit  solution  of  y  as  a  function  of  z,  we  expand  cosh  in  Taylor 
series  around  zero,  and  keep  only  the  first  two  terms  of  the  expansion.  Then, 
Eq.  (1.6)  can  be  rewritten  as 
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(1.7) 


The  results  of  the  integral  can  be  put  into  the  following  form 

Z  cos  K  z  (K  c)'1  sin  K  z  y 

y  1  .  f  °  O  o  i  r70  i 

Z  ‘  l-K  c  sin  K  z  cos  K  zJ  >-v  J 

V  o  o  o  yo 
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where  Kq  «  $  t*'C  wavenuraber  tlie  betatron  oscillations.  4n  equally 

convenient  form  for  y  is 


y  -  yB  cos  (KqZ  -  ♦„), 


where 


(2( 


Vyo/ci2  £  "2 i 1/2 

P  ,k 

ol  w 


is  the  amplitude  of  the  betatron  oscillation  and 


c‘»s'1(yo/yB) 


(1.9) 


(I. 10) 


(I. 11) 


is  the  initial  phase  of  the  betatron  oscillation 


